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ABSTRACT 



> 

\q We present equivalent widths of the [OII]A3727 and Ha nebular emission lines 

for 77 brightest cluster galaxies (BCGs) selected from the 160 Square Degree ROSAT 
X-ray survey. We find no [OII]A3727 or Ha emission stronger than —15 A or —5 A, 
respectively, in any BCG. The corresponding emission line luminosities lie below L ~ 
6 x 10 40 erg s _1 , which is a factor of 30 below that of NGC1275 in the Perseus cluster. A 
comparison to the detection frequency of nebular emission in BCGs at z ^ 0.35 drawn 



from the Brightest Cluster Survey (Crawford et al. 1999) indicates that we should 



have detected roughly one dozen emission-line galaxies, assuming the two surveys are 
selecting similar clusters in the X-ray luminosity range 10 42 erg s _1 to 10 45 erg s _1 . 
The absence of luminous nebular emission (ie., Perseus-like systems) in our sample is 
consistent with an increase in the number density of strong cooling flow (cooling core) 
clusters between z = 0.5 and today. The decline in their numbers at higher redshift 
could be due to cluster mergers and AGN heating. 



Subject headings: galaxies: clusters: cooling flows: intergalactic medium; X-rays: galax- 
ies: clusters 
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INTRODUCTION 



Galaxy clusters are the largest virialized arrangements of dark matter and baryons in the 
Universe. Most of the baryons in clusters reside not in galaxies, but in hot, diffuse atmospheres 
that are easily detected in the X-ray band. X-ray atmospheres are confined by deep gravitational 
wells, so catalogs of X-ray-selected clusters are nearly insensitive to false concentrations of galaxies 
projected along the line of sight. Furthermore, X-ray luminosities provide a reliable estimate of 
cluster mass. For these and other reasons, X-ray cluster surveys are ideally suited for selecting 
and studying the systematic properties of clusters and their galaxy populations (Rosati, Borgani, 
& Norman 2002). 

Brightest cluster galaxies (BCGs) reside near the X-ray centroid where a cooling flow (cool 
core) is often found (Fabian 1994). A cooling flow is characterized by a bright, central cusp of 
X-ray emission from gas with a radiative cooling time that is shorter than 1 Gyr. BCGs centered in 
strong cooling flows usually harbor bright nebular emission, star formation, and cold molecular gas 
(Heckman et al. 1989; Edge 2001; Salome & Combes 2003; Donahue et al. 2000). The connection 
between cold gas, star formation, and cooling flows has been strengthened by the recent discovery 
of a star formation threshold in BCGs. High star formation rates and bright nebular emission are 
found preferentially in BCGs centered in atmospheres with central cooling times shorter than ~ 0.5 
Gyr (Rafferty et al. 2008, Cavagnolo et al. 2008, Hudson et al. 2010). Therefore, the presence or 
absence of nebular emission in BCGs signals the presence or absence of an active cluster cooling 
flow (see also Edwards et al. 2007 and Balogh et al. 2002). 

We searched for nebular emission in a sample of 77 BCGs drawn from the 160 Square Degree 



(160SD) ROSAT X-ray cluster survey ( |Mullis et aLl|2003[ |Vikhlinin et al.||1998[ ). We present 
nebular emission-line strengths of [OII]A3727 and Ha, and we compare them to BCGs from the 



Brightest Cluster Survey (BCS) (Ebeling et al. 1998 Crawford et al. 1999). Absorption lines of 
H/3 and Mglb , and the molecular band strength of Mg2 are presented. No nebular emission or 
evidence of star formation is detected in any of our BCGs, which has interesting implications for 
the evolution of galaxy clusters and cooling flows. 



We assume a concordance ACDM cosmology with H D = 70 km s 1 Mpc 



2. BRIGHTEST CLUSTER GALAXY SAMPLE AND OBSERVATIONS 

2.1. The 160 Square Degree ROSAT Cluster Survey 

The 160SD survey (Vikhlinin et al. 1998; Mullis et al. 2003) yielded a catalog of 201 
galaxy clusters detected by extended X-ray emission in 647 ROSAT PSPC images. The sur- 
vey covered roughly 160 deg 2 of the sky at high fluxes, approximately 80deg 2 at a median flux 
of 1.2 x 10~ 13 erg cm -2 s _1 and 5 deg 2 at 3.7xl0 _14 ergcm _2 s _1 . The survey includes targets at 
galactic latitudes greater than 30° and with hydrogen column densities below 6 x 10 20 cm -2 . Re- 
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gions of 10° radius around the small and large Magellanic Clouds were excluded. The redshift range 
of the survey was 0.015 < z < 1.261, with a median redshift of (z) = 0.25. 



Most extended X-ray sources away from the Galactic plane are galaxy clusters (Mullis et al. 



2003) that are distinguishable from nearby elliptical galaxies and distant quasars. Sample selection 



is discussed throughly in Vikhlinin et al. (1998). All cluster candidates were imaged in the optical 
band in order to confirm their association with galaxy concentrations (McNamara et al. 2001). 
BCGs were identified by overlaying ROSAT X-ray contours onto the optical images (eg., Figure 
[I]), and spectra were taken for one or more bright galaxies nearest to the X-ray centroid of each 
cluster to measure redshifts (Mullis et al. 2003). When no obvious BCG could be identified, the 
galaxy with the brightest spectral continuum was considered the BCG. Of the 155 long slit spectra 
available from the 160SD survey (45 of the 200 cluster redshifts were found in the literature), 77 
BCG spectra were sufficiently sensitive for this study. Most of the clusters with redshifts taken from 
the literature are nearby, low luminosity systems whose redshifts were determined using absorption 
line spectra. These systems have been excluded from our study. 

The redshift range of this sub-sample is z = 0.066 — 0.7. The redshift distribution is shown 
on the left panel of Figure [2j The mean and median redshifts are z = 0.286 and (z) = 0.256, 
respectively. The X-ray luminosity distribution of our clusters lies between Lx = 0.04 — 4.27 x 
10 44 erg s -1 , as shown on the right panel of Figure[2] The sample is weighted toward low luminosity 
clusters, with mean and median values of Lx = 7.2 x 10 43 ergs -1 and (Lx) = 4.7 x 10 43 ergs -1 , 
respectively. 



2.2. Observations 



Three telescopes were used to collect spectra for this study: the University of Hawaii (UH) 
2.2 m, Keck-II 10 m and the Multiple Mirror Telescope (MMT). All 160SD survey spectroscopic 
observations include the BCG and often two or three additional galaxies. Spectra were gathered 
for 66 BCGs with the UH 2.2 m telescope using a 1.8" slit width, an effective wavelength coverage 
of 3800-9000 A, and dispersion of roughly 3.6 A pixel 1 . The spectra for six of the most distant, 
z > 0.5 BCGs were obtained with the Keck-II telescope. Keck provided an effective wavelength 
coverage of 5000-9000 A and a dispersion of roughly 2.45 A pixel -1 using a 1.5" slit width. Six BCG 
spectra were obtained with the MMT using a slit width of 1.5", a wavelength coverage between 
3600-8400 A, and a dispersion of roughly 1.96 A pixel -1 . Figure [3] shows an example spectrum of a 
BCG observed with the UH telescope that has been smoothed with a boxcar average over 5 pixels. 
The locations of spectral features measured in this study are indicated. 
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DATA REDUCTION AND ANALYSIS 



The data reduction was performed using IRAF. The spectra were de-biased, flat-fielded, back- 



grounds were subtracted, and the wavelength scale was measured. Details are given in Vikhlinin et 



al. 


(1998 


) and 


Mullis et al. 


(2003 



(2003). The spectroscopic redshifts were determined using the 4000 A 
break and absorption features of Ca II H (3933.68 A) and K (3968.49 A), the G band (~4300A), 
Mg lb (~5175A), and Na Id (5889.97 A), when detected. 



3.1. Line Equivalent Widths 

Line equivalent widths (W) were measured for [OII]A3727 and Ha in emission and H/3 and 
Mglb in absorption using the formula 

W = | F ^- F ^dA, (1) 

where F c is the flux of the continuum and is the flux of the spectral line. We follow a sign 
convention where emission lines have negative equivalent widths, and absorption lines have positive 
equivalent widths. The continuum was measured using adjacent sidebands redward and blueward 
of the feature bandwidth. A linear fit was applied across the continuum sidebands. Feature band- 
widths and continuum sidebands are listed in Table[l] flO'CoimeU||l973] |McNamara & O'ConnelT 



1989t|Fabian|1992[ ). The sidebands were placed in close proximity to the feature to avoid differential 
extinction and continuum flux calibration uncertainty, and to avoid strong absorption or emission 
features. The full spectra have not been corrected for line-of-sight extinction. Local extinction 
corrections to the emission line upper limit measurements have been applied. 

Emission lines emerging from warm ionized gas are expected to be Doppler broadened due 
to their motions in the center of the galaxy. Doppler broadened [OIIJA3727 and Ha emission 



lines in cooling flow BCGs have velocity widths lying between ~ 200 — 700 km s ^Heckman et 



al. 1989). Sampling dictates that a minimum of 3 pixels be included at the line position for 
each equivalent width measurement. For our minimum dispersion of 3.42 A pixel -1 , the minimum 
achievable line sampling corresponds to 497kms _1 at Ha and 818 km s™ 1 at [OIIJA3727. At our 
maximum dispersion 1.44 A pixel -1 , velocity widths approaching 264 km s -1 for Ha and 450kms _1 
for [OII]A3727 can be resolved. Our velocity width coverage for [OII]A3727 lies between 450 — 
818kms _1 and 264 — 497 km s -1 for Ha. Our spectra are therefore reasonably well matched to 
the expected gas velocity spread in these systems. 

The estimated uncertainty in the equivalent width measurements assume that the standard 
deviation of the best linear fit to the continuum is equal to the uncertainty at the line center. The 
uncertainties given in Table 2 are la standard deviation of the best linear fit. 
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3.2. Spectral Indices 

The spectral index, I, is the ratio of the flux of the molecular band to the interpolated contin- 



uum flux in magnitudes (O'Connell 1973). We include a measurement of I for the Mg2 molecular 



band (MgH + Mg b) to compare our BCGs to local giant elliptical galaxies ( Davies et al.|1993 ). The 
molecular and continuum bandwidths for Mg2 are listed in Table[TJ To estimate the continuum, 
the average flux of each continuum sideband was used to interpolate across the molecular band- 



width (O'Connell 1973; Faber et al. 1985). The spectral index was measured using the following 



expression: 



-2.5 log 



F x (Xi 



Fx(X 2 ) 

+ (F x (\s) - F x (Xi)) 



(2) 



In this expression, F\(\2) is the mean flux of the molecular bandwidth, F\(Ai) and F\(X^) are 
the estimated mean fluxes of the continuum bandwidths, A2 is the wavelength at the bandwidth 
center, Ai and A3 are the wavelengths at the center of the continuum bandwidths ( Q 1 Connell|| 1973 ). 
The uncertainty in the spectral index measurement was estimated by applying the method of partial 
derivative error propagation and is presented in Table 3 as the la error estimate. 



4. RESULTS 



4.1. Spectral Absorption Line Indices 



Though they were selected using similar X-ray criteria, the BCGs in this sample lie, for the 
most part, at higher redshifts than those in the Crawford et al. (1999) sample of BCGs drawn 
from the ROSAT Brightest Cluster Survey. It is therefore worthwhile to compare key line indices 
for our galaxies to those for Crawford's galaxies and to nearby giant ellipticals ( Davies et al.fl993" ) 
with similar luminosities and linear aperture sizes to identify any gross differences in their stellar 
populations. In Figure[4]we present a histogram of the distribution of iMg 2 indices for our sample. 
Our sample yields a mean of iMg2 = 0.24 ± 0.07 mag. In comparison, we found the mean Mg2 
strength of the Crawford et al. (1999) sample to be 0.28 ± 0.1, which is consistent with the mean 



for our sample. The horizontal bar in Figure 4 represents the range of iM g2 from the Davies et 
al. (1993) sample of local giant elliptical galaxies (gEs). Of the BCGs plotted in Figure [4j 65% 
lie within the range of local gEs, and 97% lie within la of the mean. Our Mg2 strengths are also 
consistent with those for gEs in the Sloan Digital Sky Survey (Bernardi et al. 2003). 

The lower-left panel of Fig ure[5| s hows the distribution of Wjj/? for our sample. The range 
of Wh/3 for the Davies et al. (1993) sample is ~0.4 — 2. 8 A, while the mean of our sample is 
W H/3 = 1.57 ± 1.52 A. The BGCs fall nicely within this range. Considering the entire sample of 
BCGs for which reliable measurements of Whs are available, ~70% fall within the Davies et al 
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( 1993 ) measured range, and ~88% fall within la of the distribution. 



An issue associated with this analysis concerns the effect of metallicity gradients on our line 
strength measurements. It is well known that the mean stellar metallicity in gE galaxies and BCGs 
generally declines with radius. This implies that for a fixed slit aperture, the observed spectrum 
of a distant galaxy extends to larger radii compared to a nearby galaxy. As a consequence, the 
stellar populations in the most distant galaxies will appear to have lower mean metallicities. The 
upper-right panel of Figure [5] shows the relationship between WMgib and redshift for the sample. 
The plot of WMgib versus z shows no trend toward declining line strength in more distant galaxies. 
Therefore, we find no discernible aperture-based metallicity bias within our sample. 

The upshot of this analysis is that we detect no gross anomalies in the two absorption line 
indices studied here that would indicate the presence of unusual stellar populations in our BCGs 
with respect to our comparison sample and normal local giant elliptical galaxies. 



4.2. Lack of Nebular Emission 

Emission lines of [OII]A3727 and Ha trace cooler gas that has been ionized by hot, young 
stars and an active nucleus. The [OII]A3727 emission-line is sensitive to low levels of ionized gas 
because of its intrinsic strength and high contrast against the declining stellar UV continuum. 
However, measuring Wrorrj can be challenging given the local sensitivity to metallicity and the 
declining sensitivity of CCD detectors in the near UV dBalogh et alpOO^ . 



Because the signal-to-noise ratios of our spectra were optimized for redshift measurements 
and not for the detection of weak line emission, we shall restrict our discussions to galaxies of 
moderate to high emission line strengths, which would be easily detected. We have chosen to 
present our upper limits primarily in the form of equivalent widths, rather than line luminosity for 
practical reasons. First, equivalent width is a differential quantity that can be measured accurately 
and independently of the flux calibration. Because our spectra were taken over several years and 
with three different instruments, equivalent width measurements avoid any systematic differences 
in flux calibration between the various observing runs. Second, Crawford et al. (1999) presents 
line luminosities for the objects they detected, but they do not present upper limits. Therefore, 
a direct comparison between the two samples in luminosity would not be feasible. Furthermore, 
nebular line emission is often spatially extended, and comparisons between line luminosities would 
be hampered by unknown aperture corrections. Equivalent widths, being differential quantities 
are less affected. Crawford et al. (1995) presented equivalent widths for a large number of their 
objects, which allowed us to convert their line luminosities to equivalent widths for the purpose of 
this comparison. 

Figure [6] shows the equivalent widths of [OII]A3727 and Ha as a function of redshift. For 
objects below redshift 0.2, the [OII]A3727 feature falls in a noisy region near the edge of the 
wavelength coverage of the spectrum. For this reason, 16 Wroni measurements were excluded 
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from Figure [HJ Similarly, measurements of Wh« for nine galaxies with redshifts above 0.35 have 
been omitted from Figure [6] because they fall outside of the wavelength coverage of the spectra. 
An additional 13 measurements were excluded because of high noise or their coincidence with a 
sky line. The [OII]A3727 and Ha features are complementary tracers of nebular emission. One 
or both were measured for each BCG. 

The upper- left panel of Figure [6] shows W[qh] as a function of redshift. The mean equivalent 
width is W[on] = 0.56 ± 4.39. No emission line stronger than 3a of the mean, or ~ —15 A is 
detected. The largest [OIIJA3727 emission line strength is W [on] .max = —13.6 ± 9.8A, but this 
feature is probably spurious. The upper-right panel of Figure [6] shows Wh<* as a function of 
redshift. The increasing uncertainty in Wh« with redshift results from increased noise in the Ha 
spectral region for objects at higher redshifts. The mean for the sample is W Ha = 0.64 ± 1.76 A. 
Again, no significant Ha emission above 3cr of the mean, or ~ — 5 A is detected. In fact, most 
show Ha in absorption. 

One sigma upper limits to the line luminosities were calculated for each of the [OIIJA3727 and 
Ha emission lines. Six uncalibrated MMT spectra were excluded. Emission-line luminosities were 



corrected for line-of-sight extinction following Cardelli, Clayton, Sz Mathis (1989). The average la 

l 



line upper limits for [OIIJA3727 and Ha are L 



[Oil] 



6. 1 x 10 40 erg s" 1 and L Ha = 3.8 x 10 40 erg s" 



respectively. The median values are (L[ n]) = 3.7 x 10 40 ergs _1 and (Lna) = 2.1 x 10 40 ergs _1 . 
Table [2] lists the L[ n], Lh q , and corresponding equivalent widths, the cluster redshifts, and their 
X-ray luminosities. Table [3] lists the equivalent widths for H/3 and Mglb, as well as the Mg2 
spectral indices for each cluster. No significant [OII]A3727 emission stronger than ~ —13 A or Ha 
emission stronger than ~ —5 A has been detected in any BCG in our sample. 

The absence of strong line emission in any 160SD BCG is surprising. The strongest emission 
lines seen in cooling flow BCGs have [OII]A3727 equivalent widths that are significantly greater 
than our upper limits. For example, [OII]A3727 reaches an equivalent width of —75.9 A in the 
Abell 1795 BCG and 
measurements of < 



-39.6 A in the Abell 2052 BCG (McNamara & O'Connell 



Crawford et al. 



( 1995 ) yielded values of -58.0 A for RXJ0439.0+0520 and 



1989}. The W [0 ii] 
23. 5A 



for the A291 BCG. Objects of this nature are clearly missing from our sample. 



5. Comparison to The ROSAT Brightest Cluster Survey 



The Brightest Cluster Survey (BCS) is a compilation of 201 clusters selected from the ROSAT 
All Sky Survey (RASS), the majority of which have X-ray luminosities exceeding Lx = 10 44 ergs -1 . 



Most lie below redshift z = 0.3 ( Ebeling et al.|[T998 ). Crawford et al. (1999) searched for emission- 
lines in a sample of 177 BCGs in clusters selected from the BCS to an unabsorbed X-ray flux limit 
of 7.9 x 10~ 12 erg cm -2 s . Their sample yielded a detection fraction of 27% in Ha emission, while 
32% show emission in other lines such as [Nil] AA6548, 6584. 

The left panel of Figure [7] compares the X-ray luminosities of Crawford's BCS clusters (indi- 
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cated as "+") to our sample (indicated by "o") versus redshift. We have multiplied the 160SD 
X-ray luminosities in this figure by the factor 1.7 to reflect the passband difference between 
BCS luminosities and ours. Figure [7] shows that the highest BCS X-ray luminosities at roughly 
Lx = 4 x 10 45 ergs -1 exceed the most luminous 160SD cluster by about five times. Nevertheless, 
there is substantial overlap in the X-ray luminosities of clusters from both surverys. Apart from 
a few low luminosity clusters, we find that at a fixed X-ray luminosity, the 160SD clusters lie at 
higher redshifts compared to their BCS counterparts. 

The right panel of Figure [7] shows only the BCGs from Crawford's sample with line emission at 
levels above our detection threshold. Crawford et al. found a roughly constant detection frequency 



of line emitting BCGs with redshift and luminosity (Crawford et al. 1999). Accordingly, this 



diagram shows many line emitting BCGs in clusters with comparable X-ray luminosities to 160SD 
clusters. Therefore, the absence of line emission in 160SD BCGs is not entirely due to their 
association with lower luminosity clusters. One distinguishing characteristic is that BCS clusters of 
a given X-ray luminosity lie at lower redshifts compared to their 160SD counterparts. We explore 
this below. 

The histograms in Figure [8] show the X-ray luminosity distributions of the BCS and 160SD 
clusters discussed here. BCS clusters harboring emission line BCGs are hatched. The 160SD X- 
ray luminosities have been corrected to the RASS bandpass. The [OII]A3727 and Ha detections 
from the BCS are shown in the left and right columns, respectively. The top row summarizes 
all detections. The second row shows detections stronger than ~ — 10 A for [OII]A3727 and Ha 
respectively, while the bottom row shows detections stronger than ~ —20 A. 

In order to make this comparison, we converted Crawford's emission line luminosities to equiv- 



alent widths. The conversion was made using several BCGs with line fluxes from Crawford et al 



(1999) found in common with Wjqii] measurements from McNamara Sz O'Connell (1989). The 
average ratio of equivalent width to flux was multiplied by the [OII]A3727 line flux reported by 



Crawford et al. (1999) for each line emitting BCG. These ratios are 5.2±2.6 for [OII]A3727 and 



9.7±4.5 for Ha. Of the 177 BCS clusters presented by Crawford et al. (1999), 10 were omitted 



from Figure [8] because their X-ray luminosities were unavailable to us. Of these, four BCGs have 
Ha in emission and two have [OII]A3727 in emission. Excluding these objects, the BCS con- 
tains 44 BCGs with [OII]A3727 in emission and 49 with Ha in emission. Fifty seven percent of 
the [OII]A3727 emission-line BCGs and 67% of Ha emission- line BCGs have equivalent widths 
stronger than ~ —10 A. Thirty four percent of [OII]A3727 emission-line BCGs and 43% of Ha 
emission-line BCGs show equivalent widths stronger than ~ —20 A. 

The nebular emission detection thresholds for the BCS shown in Figure [8] lie comfortably above 
ours. Such objects, if present in our sample, would have been detected. Of the entire sample of 177 
BCS BCGs from Crawford et al. ( |1999[ ), 121 fall within the X-ray luminosity range of this sample. 
Within this luminosity range, 7% of BCS BCG have [OII]A3727 emission stronger than ~ —20 A 
and 16% have Ha emission stronger than ~ —10 A. No 160SD BCG reveals nebular emission at 
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this level. 

The right panel of Figure [7] shows the X-ray luminosities of the 160SD BCGs and BCS 
BCGs harboring nebular emission above our thresholds versus redshift. The BCGs revealing only 
[OII]A3727 emission are marked with a "O"; those with both [OII]A3727 and Ha emission are 
indicated as "+". This figure shows that at a given redshift, BCS clusters have higher average 
X-ray luminosities. 

About 30% of BCS BCGs with cluster X-ray luminosities greater than Lx = 10 45 erg s _1 have 
emission lines. The fraction falls to 18% within the X-ray luminosity range 10 44 ergs _1 > Lx > 
10 45 ergs -1 . Similarly, 16% of clusters with X-ray luminosities within the range 10 42 ergs _1 > 
Lx > 10 44 ergs _1 have emission lines. These figures show that while higher luminosity clusters 
show a higher frequency of emission line BCGs, the lower average X-ray luminosity of the 160SD 
sample cannot alone explain the absence of strong line emission. 

Table [4] summarizes the detection statistics for both surveys. Assuming both samples are 
equally likely to have BCGs with bright line emission within the overlapping luminosity range, 
the table shows that we should have detected line emission in roughly one dozen 160SD BCGs 
with Ha equivalent widths greater than —10 A, and with [O II] widths exceeding —20 A. The 
probability of having detected one or fewer emission line BCGs by chance is approximately one 
in ten thousand. Eight percent of BCS clusters within our X-ray luminosity range host BCGs 
with line equivalent widths stronger than —20 A. We would have expected to detect about 6 such 
systems in our sample (Figure 8, panel 3), yet we detected none. The likelihood of having detected 
one or fewer such system by chance is less than 2%. Thus, the absence of line emitting BCGs in 
the 160SD survey is statistically significant at a high level. 

5.1. Potential Systematic Errors 

A potential systematic effect that could bias our equivalent widths downward is beam dilution 
owing to larger projected linear aperture sizes with increasing redshift. Crawford et al. used a slit 
width of approximately 1.2 arcsec, which is somewhat smaller than ours. Crawford's clusters lie 
on average at lower redshifts, so that their slit width projects over a several kpc on the galaxy. In 
comparison, our slits subtended a linear width of between 6 — 9 kpc between redshifts of 0.3 and 
0.5. Furthermore, Crawford integrated along the slit out to a typical linear distance of about 10 
kpc, while we integrated along the slit to the approximate location where the galaxy brightness 
falls below the sky brightness, which is 3 — 5 times larger in radius. Because of the larger fraction 
of galaxy light in our aperture compared to Crawford's, our equivalent may be diluted with respect 
to theirs. 

We have estimated the magnitude of the effect by modeling the light profile of our galaxies with 
a Hubble law, ie., stellar surface brightness declining radially as r -2 , and integrating the stellar 
light beyond the central 10 kpc down to the sky limit. We found that the effect would bias our 
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equivalent widths downward by at most 30% — 40% with respect to Crawford's. In other words, 
it is possible that an object with Ha equivalent width of — 7A would fall below our — 5A cutoff 
because of aperture dilution. However, the detection of stronger features is not in jeopardy. In 
fact, an upper limit on the magnitude of the effect can be evaluated directly by scanning the line 
luminosity upper limits for the most distant objects shown in Table 2, which would show the largest 
effects. These upper limits lie below 10 41 erg s _1 . In comparison, the objects that are putatively 
missing from our sample would have line luminosities lying between 10 41 erg s -1 and 10 43 erg s _1 . 
They would easily have been detected. In fact, Crawford's smaller aperture will have missed flux 
from spatially extended emission line systems that we would have captured. Many emission line 
systems in cooling flows, such as Perseus, extend 30 kpc or more from the nucleus (Heckman et al. 
1989, O'Dea et al. 2010). Therefore, a dilution effect will be most pronounced in weak emission 
line galaxies with small, unresolved nebulae, which we are not concerned with in this paper. 

While the detectability of strong emission lines is not at issue, our primary concern is whether 
the X-ray selection functions in the 160SD and RASS surveys are oppositely biased toward the 
non detection and detection of cooling flow clusters in each respective survey. In a study of 10 
BCGs from the 160SD survey with X-ray luminosities below Lx = 4 x 10 43 ergs _1 [0.1-2.4 keV], 



Balogh et al. (2002) also failed to detect line emission in any BCG. They argued that the 160SD 
survey is biased against clusters with strong cooling flows because only extended sources at a faint 
X-ray flux limit were sought. This interpretation is incorrect. If the 160SD survey is biased, we 
would expect it to contain an excessive number of cooling flows. The central ~ 100 kpc or so 
of cooling flow clusters are brighter than normal clusters. They are well resolved in the 160SD 
survey, and are thus easier to detect than non-cooling flow clusters (Vikhlinin et al. 1998, Burenin 
et al. 2007). Simulations have shown that the detection probability declines precipitously only for 
clusters harboring X-ray- bright AGN near their centers, and only in the rare event that the AGN 
outshines the cluster emission (Burenin et al. 2007). Such clusters are rare at these redshifts. 

Figure 7 shows that clusters of a given X-ray luminosity in the BCS sample lie at systematically 
larger redshifts compared to their 160SD survey counterparts. If the lack of emission lines in the 
160SD survey is a physical effect and is not due to selection bias, it appears to be related to redshift. 
This interpretation implies that at a given X-ray luminosity, there are fewer well established strong 
cooling flows at higher redshift, suggesting rapid evolution between redshift one-half and today. 

Note that our relatively high detection threshold is insensitive to the existence of weak cooling 
flows in the 160SD survey. Roughly 15% of BCGs selected from the Sloan and NOAO Fundamental 
Plane cluster surveys (Edwards et al. 2007) have detectible line emission. This fraction rises to near 
unity for bona-fide cooling flow clusters. However, the Ha equivalent widths of most are weaker 
than our — 5A threshold. Nevertheless, 28 of their 328 BCGs (8%) drawn from the Sloan survey 
have Ha equivalent widths stronger than — 5A, and all of them lie at redshifts at or below ~ 0.1. A 
direct comparison between our survey and the Edwards et al. (2007) survey is infeasible because of 
the different selection methods. However, it is clear that were similar BCGs present in our survey, 
we would have detected a substantial fraction of them. The absence of any BCG with strong line 
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emission in our sample bolsters the case for an evolutionary effect. 

6. DISCUSSION 

We searched for nebular emission in a sample of 77 BCGs from the 160SD X-ray cluster 
survey. No [OII]A3727 or Ha emission stronger than —15 A and —5 A, respectively, was detected 
in any BCG. A comparison between our sample and the BCS clusters of similar X-ray luminosity 
suggests that we should have detected roughly one dozen emission-line BCGs stronger than our 
detection threshold. 160SD clusters of a given X-ray luminosity lie at systematically higher redshifts 
compared to BCS clusters of the same luminosity. This may indicate an increase in the number 
density of strong cooling flows between z = 0.5 and today. Because of our relatively high detection 
thresholds, our result does not imply that cooling flows don't exist at the redshifts surveyed in the 
160SD survey (see Bauer et al. 2005). However, we find a significant decrease in the number of 
strong cooling flows like the Perseus Cluster. Similar conclusions were reached using independent 
methods by Santos et al. (2008), Vikhlinin et al. (2007) and by Voevodkin et al. (2010), who also 
found fewer cooling flow clusters at higher redshifts. 

The decline may be related to the destruction of cooling flows by mergers, which occurred 
more frequently than they do today (Santos et al. 2008). Weaker mergers may also disturb X-ray 
atmospheres and miscenter the BCG, which could prevent a strong cooling flow from forming. AGN 
feedback is thought to be responsible for quenching cooling flows in nearby clusters (Peterson & 
Fabian 2006; McNamara & Nulsen 2007), but this is less clear in higher redshift clusters lacking 
the signatures of cooling flows. However, in a study of radio AGN located in the 400 Square Degree 
X-ray cluster survey (Burenin et al. 2007), Ma et al. (2011, in preparation) found that the power 
injected by radio sources in galaxies lying within 250 kpc of cluster centers is significant compared 
to the power radiated by the hot atmospheres. The selection function of the 400 Square Degree 
Survey is nearly identical to the 160 SD survey, but it covers a much larger area of the sky. The 
Ma et al. study suggests that AGN are able to deposit enough energy in these clusters to delay or 
prevent a massive cooling flow from forming at redshifts of a few tenths or so. 
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Fig. 1 — An optical photometric image (10' x 10') of cluster RXJ0056. 9-2740, with ROSAT PSPC 
X-ray contours overlaid. The cluster is at a redshift of z = 0.116 with an X-ray luminosity of 
Lx = 1.6 x 10 44 ergs _1 . The bright optical source at the top of the image is a foreground star. 
The brightest galaxy located at the X-ray centroid of the cluster is clearly discernible from the 
surrounding cluster members. 
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Fig. 2. — (Left) The redshift distribution of the 77 BCGs represented in this sample. The sample 
ranges in redshift from 0.066 — 0.7, with a median redshift of (z) = 0.256 and a mean redshift of 
z = 0.286. (Right) The X-ray luminosity distribution on a logarithmic scale of the host clusters 
for the sample. The X-ray luminosity for the clusters range from 0.04 — 4.27x 10 44 erg s _1 , with a 
mean X-ray luminosity of Lx = 7.2 x 10 43 ergs -1 and a median of (Lx) = 4.7 x 10 43 ergs -1 . 
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Fig. 3.— A longslit spectrum of the BCG from cluster RXJ1343.4+5547 at z = 0.069 with L x = 
4.0 x 10 42 ergs _1 , taken with the University of Hawaii 2.2 m telescope. The vertical dashed lines 
mark the location of emission and absorption features, with one large atmospheric O2 absorption 
feature. The spectrum has been smoothed using a boxcar average over 5 pixels 
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Fig. 4. — A histogram of spectral indices for the Mg2 molecular band for the sample of BCGs. 
The boldface bar plotted above the histogram denotes the range of Mg2 values indicative of gE 

. The mean of the sample (Ijvig2 = 0.24 ± 0.07 mag) falls central to the 
the range of local gEs. 



galaxies (Davies et al. 1993 
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Fig. 5. — (Upper- left) The equivalent width of H/3 versus redshift for the sample of BCGs. 
Compared with local gEs, the BCGs show similar H/3 equivalent widths. (Lower-left) The dis- 
tribution of H/3 equivalent widths for the BCGs, centered about a mean equivalent width of 
Wh/3 = 1.57 ± 1.52 A. (Upper-right) The equivalent width of Mglb versus redshift, which shows 
no significant gradient with redshift. (Lower-right) The distribution of measured Mglb equivalent 
widths versus redshift. The adopted convention is that negative values indicate emission, while 
positive values indicate absorption. 
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Fig. 6. — (Upper-left) The equivalent width of [OII]A3727 versus redshift for the BCGs of this 
sample. The mean of [OII]A3727 equivalent width is Wron] = 0.56 ± 4.39, and no significant 
[OII]A3727 emission stronger than ~ —15 A was detected in the sample. (Lower-left) The distri- 
bution of [OII]A3727 equivalent widths among the BCGs. (Upper-right) The equivalent widths 
of Ha versus redshift, where the mean is W Ha = 0.64 ± 1.76 A, and no significant Ha emis- 
sion stronger than ~ — 5A has been detected in the sample. (Lower-right) The distribution of 
Ha equivalent widths among the BCGs. The adopted convention is that negative values indicate 
emission, while positive values indicate absorption. 
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Fig. 7. — (Left panel) A comparison of the BCS sample of Crawford et al. (1999), marked with 
'+', and this sample, marked with 'o', with X-ray luminosity plotted as a function of redshift for 
the BCGs. (Right panel) The emission-line BCGs of the BCS sample plotted with cluster X-ray 
luminosity as a function of redshift. The BCGs that yield Ha emission stronger than —10 A are 
marked with a 'O', while BCGs that yield [OII]A3727 emission stronger than —20 A are filled with 
a '+'. The BCGs of this sample are also plotted and marked with 'o'. 
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Table 1. The feature band widths and continuum sidebands used to measure line strengths. 



Feature 


Line Center (A) 


Bandwidth (A) 


Continuum Sidebands (A) 


[Oil] 


3727 


3721.00-3732.00 


3692.00-3722.00, 3732.00-3762.00 


Ha 


6562 


6557.00-6568.00 


6527.00-6557.00, 6567.00-6597.00 


R/3 


4861 


4849.50-4877.00 


4829.50-4848.25, 4878.25-4892.00 


Mgib 


5175 


5162.00-5193.25 


5144.50-5162.00, 5193.25-5207.00 


Mg 2 t 




5156.00-5197.25 


4897.00-4958.25, 5303.00-5366.75 



f Mg2 is a molecular bandwidth measured using a spectral index in magnitudes. 
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Fig. 8. — A set of histograms that show the X-ray luminosity distribution of BCGs for the Brightest 
Cluster Sample (BCS), this sample (labeled as 160SD), and the BCS emission line galaxies. For 
comparison purposes, the X-ray luminosities of this sample have been corrected to the RASS band- 
width [0.1-2.4 keV]. The left and right columns are designated for BCGs that show [OII]A3727 and 
Ha emission respectively. The second row shows detections stronger than ~ — 10 A for [OII]A3727 
and Ha , while the third row shows detections stronger than ~ —20 A. 
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Table 3. The absorption line strengths for the BCGs. 



Object 


W H/3 


WMglb 


lMg 2 




(A) 


(A) 


(mag) 



RXJ0041. 1-2339 


2.24 


± 


0.64 


3.30 


± 


0.94 


0. 


,26 


± 


0.11 


RXJ0050.9-0929 


1.07 


± 


0.42 


2.97 


± 


0.73 


0. 


.31 


± 


0.07 


RXJ0056. 9-2213 


2.32 


± 


0.90 


3.79 


± 


0.72 


0. 


.29 


± 


0.13 


Rxjono.s+igss* 


1.90 


± 


0.85 








0. 


.27 


± 


0.10 


RXJ0122.5-2832 


2.32 


± 


0.41 


3.95 


± 


0.60 


0. 


.29 


± 


0.11 


RXJ0124.5+0400 


2.17 


± 


0.92 


4.19 


± 


1.03 


0. 


.27 


± 


0.14 


RXJ0142.8+2025 


2.24 


± 


0.46 


3.09 


± 


0.91 


0. 


.24 


± 


0.12 


RXJ0144.4+0212 


0.61 


± 


0.47 


-0.54 


± 


0.78 


0. 


.21 


± 


0.06 


RXJ0159.3+0030 


1.63 


± 


1.00 


0.05 


± 


0.69 


0. 


.14 


± 


0.08 


RXJ0206.3+1511 


1.11 


± 


1.02 


3.86 


± 


1.21 


0. 


.25 


± 


0.16 


RXJ0206.8-1309* 


0.51 


± 


1.43 








0. 


.28 


± 


0.11 


RXJ0258.7+0012 


2.32 


± 


1.33 


1.93 


± 


1.30 


0. 


.25 


± 


0.16 


RXJ0259.5+0013 


2.38 


± 


0.76 


5.00 


± 


1.11 


0. 


.25 


± 


0.14 


RXJ0351.6-3649 


2.29 


± 


0.90 


4.78 


± 


0.99 


0. 


.22 


± 


0.18 


RXJ0506.0-2840 


2.37 


± 


1.13 


0.13 


± 


1.79 


0. 


.22 


± 


0.12 


RXJ0521. 1-2530 


1.13 


± 


1.77 


5.99 


± 


3.26 


0. 


.38 


± 


0.41 


RXJ0522.2-3625 J 


1.63 


± 


1.34 








0. 


60 


± 


0.60 


RXJ0826.4+3125 


2.54 


± 


0.99 


3.25 


± 


0.75 


0. 


.24 


± 


0.11 


RXJ0841. 1+6422 


-0.99 


± 


1.64 


5.74 


± 


2.51 


0. 


.25 


± 


0.28 


RXJ0842.8+5023 


1.98 


± 


1.75 


3.20 


± 


2.04 


0. 


.20 


± 


0.29 


RXJ0852. 5+1618 


1.57 


± 


1.38 


4.56 


± 


1.41 


0. 


.34 


± 


0.24 


RXJ0858.4+1357 


1.68 


± 


0.45 


2.55 


± 


0.79 


0. 


19 


± 


0.11 


RXJ0907.2+3330 


2.41 


± 


0.92 


1.90 


± 


1.56 





.24 


± 


0.25 


RXJ0921. 2+4528 


0.82 


± 


1.22 


1.35 


± 


0.84 


0. 


16 


± 


0.19 


RXJ0926.6+1242* 


3.71 


± 


2.61 








0, 


.18 




0.51 


RXJ0943. 7+1644 


2.59 


± 


0.76 


4.97 


± 


1.58 


0. 


.31 


± 


0.25 


RXJ0958.2+5516 


-0.33 


± 


1.84 


7.49 


± 


1.70 


0. 


.34 


± 


0.28 


RXJ1013.6+4933 


1.49 


± 


1.38 


3.99 


± 


1.03 


0. 


32 


± 


0.13 


RXJ1015. 1+4931 


1.77 


± 


1.48 


-0.15 


± 


0.95 


0. 


.22 


± 


0.17 


RXJ1036.1+5713 


6.68 


± 


1.44 


1.55 


± 


1.31 


0. 


.08 


± 


0.13 


RXJ1049. 0+5424 


2.64 


± 


0.36 


1.59 


± 


0.67 


0. 


16 


± 


0.10 


RXJ1117.2+1744 


2.24 


± 


1.43 


1.42 


± 


1.66 


0. 


03 


± 


0.12 


RXJ1117.5+1744t 








2.88 


± 


0.68 


0. 


.24 


± 


0.16 


RXJ1120.9+2326 t 








1.17 


± 


0.77 


0. 


.28 


± 


0.14 


RXJ1123.1+1409 


2.09 


± 


0.93 


1.38 


± 


1.58 


0. 


23 


± 


0.14 


RXJ1124.0-1700 


-0.32 


± 


1.65 


0.43 


± 


1.24 


0. 


17 


± 


0.23 


RXJ1124.6+4155 


3.88 


± 


1.19 


2.72 


± 


0.92 


0. 


.20 


± 


0.11 



Table 3 — Continued 



Object 


w H/3 
(A) 


WMglb 

(A) 


lM g2 

(mag) 


i ~> v ti 1 or n i 01 01 

RXJ 1135.9+2131 


n ia 


4 


n 54 

u. o^± 


9 67 4- 

i.UI 31 


n 48 


U.24 ± 


0.11 


Ti V T1 1 A O n i O 1 A A 

RXJ 1142.0+2144 


zl 91 


4 


1 ^n 

1 . OU 


1 R1 4 

O.U1 31 


1 ^1 

1.01 


U.41 ± 


0.16 


Ti "V T1 1 AH A i OOP" ,1 

RXJ 1146.4+2854 


n sn 

u.ou 


4 

33 


n 7Q 

U. / a 


1 9^ 4 

l.ZO 31 


n 49 


U.23 ± 


0.14 


T> V T1 1 CO 1 i r roi 

RXJ 1158.1+5521 


1 ^fi 


4 


n sn 

u.ou 


9 zL8 4 

Z.^O 31 


u.uo 


U.26 ± 


0.10 


T > -\7" T 1 nn n non^ 

KXJ12U0.9-U.327 


n ^8 


4 


1 11 

1.00 


9 ^ 4 

Z.OO 31 


1 71 


U.22 ± 


0.17 


T > -\7" T 1 c\c\c r n ^7 1 A * 

KAJ 1206.5-0744 












U.U5 ± 


0.05 


t> v ti 01 r i no r 

KAJ 1213.5+0253 


4 97 


4 

31 


1 


1 n^ 4 

O.UO 31 


1 87 
1.0 1 


U.19 ± 


0.16 


T> V T1 01 O A i 

RXJ 1218.4+3011 


9 

z.uo 


4 

31 


1 ns 

l.UO 


A 71 4 


1 fin 

l.UU 


0.22 ±0.21 


t> v ti 001 /i i a r\~\ a 

KAJ 1221.4+4918 


1 8 
u. ±0 


4 

3Z 


n 44 


6 87 4 

U.0 1 31 


n qq 

U. UCJ 


0.13 ±0.28 


TiVTIOOT £? 1 O/'OO 

KAJ 1237. o+2b«i2 




4 

31 


n 89 

U.OZ 


Q CI X 
O.Ol 31 


1 99 


0.16 


± 


0.11 


T> V T10C/1 i O C A C 

KAJ 1254.6+2545 


n 98 


4- 

3Z 


1 1 1 

1 . lO 


4 47 4 

1. 4: 1 3C 


1 11 
l.OO 


0.23 


± 


0.25 


KAJ 1254.8+2550 


1 RQ 


4 

31 


n 8Q 


1 9fi 4 

O.ZU 31 


n 71 

u. / 


0.24 


± 


0.08 


Ti v ti or/ 1 n 1 r r 

KXJ1256. 0+2556 


1 nn 


4 

31 


U.Uc/ 


^ ztn 4 

O.^U 31 


n 8zt 


0.24 


± 


0.13 


T> v ti om t i 1 nrn 

KAJ 1301.7+1059 


1 88 
1.00 


4- 

31 


z. 10 


zt ^n 4 

^.JU 31 


1 8Q 
i.oy 


0.30 


± 


0.16 


t>vti on n i 0000 

KAJ 1309.9+3222 


9 71 


4 

31 


1 zt 1 ^ 

1 .^lO 


zt fizt 4 

31 


9 9zl 


0.21 


± 


0.36 


t> v ti 007 i oor /I 

KAJ 1337.8+3854 


1 R8 
-i.uo 


4 

31 


1 ^fi 

l.OU 


n fin 4 

U.UU 31 


1 n8 

l.UO 


0.24 


± 


0.11 


T> V T1 O/IO O i a n 

KAJ 1342.8+4028 


n ^n 

u.ou 


4 

31 


U. c/O 


1 fiQ 4 
o.uy 31 


1 79 


0.28 


± 


0.36 


T> V T1 O/IO A i rr r ^ T 

KAJ 1343.4+5547 


n QD 

U. c/U 


4 

31 


n ^8 

u.oo 


1 9Q 4 
o.zy 31 


n 81 

U.Ol 


0.27 


± 


0.10 


T> "V T1 OCT /( O i pni T 

KAJ 1354.8+6917 


1 1 ^ 
-J — LO 


4- 

31 


U. 1 


^ 7^ 4 

O./O 3 


1 

1.00 


0.40 


± 


0.17 


T> V T 1 ^OO C\ 1 /"'/lOO 

KAJ 1438.9+6423 


1 "39 
o.oz 


4 

3Z 


n 6Q 


9 ^n 4 

i.OU 31 


n on 


0.25 


± 


0.15 


Ti v ti r n n n i 0000 

KXJ1500. 0+2233 


1 

1 . c/O 


4 

31 


n Q7 

u. y / 


1 99 4 

l.ZZ 31 


1 9n 

l.ZU 


0.15 


± 


0.11 


Ti "V T 1 r 1 r r 1 1 a c\ 

KXJ1515.5+4346 1 


9 49 


4 

31 


1 8^ 
1.00 


O.UO 31 


1 4^ 
1.^0 


0.25 


± 


0.10 


Ti v ti n r c* i /lorn 

KXJ 1515.6+4350 


1 riR 

l.UO 


4 

31 


n ^8 

u.oo 


9 Qfi 4 

Z. jU 31 


n 74 

u. / ^± 


0.21 


± 


0.10 


Ti v ti ro/i i nnr^ 

KXJ 1524.6+0957 


9 18 

L. lO 


4 

31 


n fi/L 

U.U^i 


q no + 

o.uy 31 


1 1 7 

J — L / 


0.17 


± 


0.17 


KXJ 1537.7+1200 


1 8zt 
1.0^ 


4 

31 


n aq 


9 n7 4 

Z.U / 31 


n 79 

U. / L 


0.25 


± 


0.06 


Ti V t 1 r /i n i 1 /| /I r 

KXJ 1540.8+1445 


1.61 


± 


1.96 


6.U7 ± 


2.14 


0.10 


± 


0.38 


Ti v t 1 r /i t i onr/ 1 

KXJ 1547.3+2056 


0.79 


± 


U.72 


1.01 ± 


1.15 


0.26 


± 


0.12 


Ti v ti rro i on 1 

KXJ 1552.2+2013 


1.88 


± 


1.18 


5.7U ± 


1.22 


0.32 


± 


0.20 


RXJ1630. 2+2434 


1.41 


± 


U.86 


4.42 ± 


U.49 


0.38 


± 


0.17 


RXJ1642.6+3935 


-0.41 


± 


2.74 


-U.9U ± 


1.71 


0.03 


± 


0.29 


RXJ1659.7+341U 


1.34 


± 


U.86 


4.24 ± 


U.88 


0.27 


± 


0.16 


RXJ1722.8+41U5 


-0.26 


± 


1.36 


2.47 ± 


1.9U 


0.32 


± 


0.22 


RXJ1729. 0+7440 


5.10 


± 


2.UU 


-U.44 ± 


2.57 


0.23 


± 


0.36 


RXJ 1746. 4+6848 


U.U4 


± 


2.U3 


2.22 ± 


1.84 


0.23 


± 


0.25 


RXJ2202.7-1902 


2.19 


± 


1.11 


3.43 ± 


1.24 


0.22 


± 


0.20 


RXJ2212. 6-1713 


1.45 


± 


U.57 


3.U6 ± 


U.44 


0.27 


± 


0.09 


RXJ2213. 5-1656 


2.15 


± 


1.35 


4.26 ± 


1.U6 


0.26 


± 


0.15 
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Table 3 — Continued 



Object 


w H/3 


WMglb 


lMg 2 




(A) 


(A) 


(mag) 


RXJ2257.8+2056 


2.62 ± 0.66 


4.04 ± 1.04 


0.24 ± 0.10 


RXJ2328.8+1453 


-1.44 ± 1.15 


3.94 ± 2.30 


0.30 ± 0.24 


RXJ2348.8-3117 


1.05 ± 0.38 


3.15 ± 0.92 


0.23 ± 0.14 



t H/3 line is coincident with an O2 sky absorption line 

* Mglb is coincident with an O2 sky absorption line 

* Mglb spectrum show local disturbance 



Table 4. The number of BCGs per X-ray luminosity bin for the BCS and this sample, as well as 

the frequency of detection for BCGs in the BCS. 



Lx Range 


160SD 


BCS 


BCS 


Fraction Detected 


(ergs" 1 ) 


BCGs 


BCGs 


Emission-line 


in BCS 








BCGst 




L x > 10 45 





33 


10 


30% 


10 44 < L x < 10 45 


29 


103 


18 


18% 


10 42 < L x < 10 44 


48 


31 


5 


16% 



^Emission- line BCGs with equivalent width estimates greater than ~ 
-20A for [OII]A3727 and ~ -IOA for Ha. 



